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ABSTRACT 

The relevance of the standing accretion shock instability (SASI) compared to neutrino-driven con- 
vection in three-dimensional (3D) supernova-core environments is still highly controversial. Studying 
a 27 M© progenitor, we demonstrate, for the first time, that violent SASI activity can develop in 
3D simulations with detailed neutrino transport despite the presence of convection. This result was 
obtained with the Prometheus- Vertex code with the same sophisticated neutrino treatment so 
far used only in ID and 2D models. While buoyant plumes initially determine the nonradial mass 
motions in the postshock layer, bipolar shock sloshing with growing amplitude sets in during a phase 
of shock retraction and turns into a violent spiral mode whose growth is only quenched when the infall 
of the Si/SiO interface leads to strong shock expansion in response to a dramatic decrease of the mass 
accretion rate. In the phase of large-amplitude SASI sloshing and spiral motions, the postshock layer 
exhibits nonradial deformation dominated by the lowest-order spherical harmonics (i = 1, m = 0, ±1) 
in distinct contrast to the higher multipole structures associated with neutrino-driven convection. We 
find that the SASI amplitudes, shock asymmetry, and nonradial kinetic energy in 3D can exceed those 
of the corresponding 2D case during extended periods of the evolution. We also perform parametrized 
3D simulations of a 25 Mq progenitor, using a simplified, gray neutrino transport scheme, an axis-free 
Yin- Yang grid, and different amplitudes of random seed perturbations. They confirm the importance 
of the SASI for another progenitor, its independence of the choice of spherical grid, and its preferred 
growth for fast accretion flows connected to small shock radii and compact proto-neutron stars as 
previously found in 2D setups. 

Subject headings: supernovae: general — hydrodynamics — instabilities — neutrinos 
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1. INTRODUCTION 

Today there is a universal consensus in supernova the- 
ory that the neutrino-driven explosion mechanism hinges 
vitally on the supporting action of multidimensional hy- 
drodynamic instabilities (with the notabl e exception of 
low-m ass progenitors with O-Ne-Mg core; iKitaura et al.l 
2006). Two different hydrodynamical instabilities can 
potentially operate in the accretion flow between the 
stalled supernova shock and the forming neutron star: 
On the one hand, neutrino heating establishes a neg- 
ative entropy gradient in the gain layer, which is un- 
stable to buoyancy-dr iven convect ion, a fact that has 
long been recognized (Bcthg I1990D and was confirmed 
already b y the first generation of multi-dimensional sim- 
ulations dHerant et all 119921: IBurrows fc Frvxell 119921: 



Herant et alj|1994t IBurrows et alj|1995t Uanka fc Mulled 
19961: iMuller fe Jankal 11997(1 . Later, purely hydrody- 



namical simulations of adiabatic accretio n flows in the 
down stream region of a stalled shock (Blondin et al. 
2003) revealed the existence of a distinctly different 
phenomenon, the "standing-accretion shock instability" 
(SASI), which involves large-scale I = 1 and £ = 
2 sloshing motions of the shock front and in three 
dimensions can also develop £ = l.ro = ±1 spi- 
ral modes (iBlondin fc Mezzacappal 120071 : llwakami et al.1 
120091 IFernandezI |2010( 1~ The characteristic oscillatory 
growth of this instability is thought to be mediated by 
an amplifying advective-acoustic cycle ( Foglizzo 120021 : 
iFoglizzo et al][200l 120071 IGuilet fc Foglizzol 12012 see- 
however. IBlondin fc Mezzaca ppa 20061 and lLaminal2007l 



for the interpretation as a purely acoustic cycle) . The hy- 
pothesis of an advective-acoustic process receives strong 
support by a detailed analysis of the mode frequencies 
(jScheck et all [20081: IGuilet fc Fogli^ol [20111 and by a 
shallow water analogue, the SWASI experiment, in which 
the water flow between a circular reservoir and a central 
tube as sink mimics t he accretion flow th at feeds the 
nascent neutron star ([Foglizzo et al.|[2012l ). With the 
role of the shock and of acoustic waves being played by 
a hydraulic jump and surface water waves, respectively, 
similar sloshing and spiral motions are observed as in the 
hydrodynamical simulations of collapsing stellar cores. 

Violent hydrodynamic mass motions in the postshock 
layer do not only have the potential to improve the 
heating conditions in the supernova core by pushing the 
shock farther out and by prolonging the dwell time of 
the accreted matter in the heating region, which has 
supportive consequences for the supernova explosion 
dHerant et al.lll994tlBurrows et al]|1995UJanka fc Mullerl 



1996; 



2009; 
[2010b : 



Fryer fc Hegerl 120001: iMarek fc Jankal 
I Murphv fc Burrows! 120081: INordhaus et al.l 
lHanke et all 120121: ISuwa et al.l I2010H . Convection 
and the SASI also create seed asymmetries that de- 
termine !_the_eje£ta_mcaTjffology_during the explosion 
phase (|Kifonidis et al.l 120031 l20ull . Moreover, large- 
scale anisotropics in the ejecta may yield a natural 
explanation for the observed h i gh ki c k veloci t ies of 
young pulsars ( Janka fc Mullerl 119941: IHerantl 1 1995 : 



Burrows fc Haves! 119961: IScheck et al l 12004 12006 



Wongwathanarat et al.1 l2010bl: INordhaus et all [2010a, 



2 



Hanke et al. 



|20H iWongwathanarat et all [20121 ) . The I = 1, \m\ = 1 
spiral mode of the SASI could also provide an effec- 
tive me ans for spinning up t h e pro to-neutron star 
(PNS) ([Blondin fc Mezzacappal l2007f l. Moreover, 
large-amplitude SASI motions were found to stir strong 
g-mode activity in the neutron st ar surface, which 
leads to gravitation al-wave emission (|Marek et alJ :2009: 
I Murphy et alj 120091 ) . Quasi-periodic variations of the 
accretion flow onto the neutron star due to SASI shock 
oscillations may also cause fluctuations of the radiated 
neutr i no luminosities and mean energies (jMarek et al.l 
[20(51 IQtt et all [200l IBrandt et all 12011 . The cor- 
responding neutrino-signal fluctuations on millisecond 
timescales could be measured for a galactic supernova 
with good time resolution as, e.g., pro vided by the Ice 



Cube detector (jLund et al.ll2010L l2012t ) and could serve 
as val uable probe of neutrino properties (Ell is et al.1 
l2012allbT ) as well as of the dynamics in the supernova 
core. 

Up to now most investigations of convection and the 
SASI in core-collapse supernovae have relied on axisym- 
metric (2D) simulations. Three-dimensional (3D) mod- 
els based on various approximations for treating neu- 
trino heating and cooling in the supernova core have 
only recently become available, but have already sparked 
a controversy about the development a nd mutual in- 
teract ion o f the two instabilitie s in 3D. iBurrows et al.1 
(2012i) and IDolence et all (|2013D . who conducted simu- 
lations using a simple light-bulb neutrino scheme, were 
rather outspoken in classifying the SASI as a subdom- 
inant phenomenon in the presence of neutrino heat- 
ing. They argued that the violent sloshing motions 
seen in 2D neutrino-hydrodynamics simulations may be 
an artifact of the artificial symmetry assumption and 
are actually n othin g but 2D convection in disguise. 
iMurphv et al.l ([20121 ) also noted that nonlinear convec- 
tion theory appears to explain the 3D flow properties 
of their models without the need of invoking the SASI 
as an additional instability. At first glance, these find- 
ings seem to be in line with other 3 D studies rely- 
ing on a similar light-bulb methodology ( Iwakami et al.| 
20081 : iNordhaus et all l2010bl: lHanke et al.l l2012t iCouch 



20121 ) or on a gray neutrino transport approximation 
with chosen neutr ino luminosities imposed at an in- 
ner grid boundary (IWongwathanarat et al.ll2010bl. [20121 : 
IMuller et al~ll2012cD . In contrast. iMiiller et al.l (|2012al ) 
demonstrated by self-consistent, 2D, general relativis- 
ts (GR) simulations with sophisticated transport that 
genuine SASI activity remains possible for sufficiently 
small shock stagnation radius (caused by high mass ac- 
cr etion rates in the par ticular case of a 27 Mq progenitor 
of iWooslev et al.ll2002t ) . This suggests that details of the 
conditions may decide about the growth of the SASI, 
and that these conditions may not only depend on the 
properties of the progenitor star but also on the exact 
behavior of the stalled shock, which again depends on 
a reliable treatment of the neutrino physics. The men- 
tioned 3D models might simply have missed the sweet 

spot for SASI growth in parameter space. 

Rec ent Newtonian simula tions by [T akiwak j et al.l 
(I2012D an d GR s imulations bv lKuroda et al.1 ([20121 ) arid 
IQtt et all (|2012f ) are first, tentative steps to higher so- 
phistication in 3D models, but these works were still 
focused on a few progenitors only and employed crude 



neutrino transport methods with various simplifications 
concerning the description of neutrino propagation, of 
neutrino intera ctions, an d of th e energy dependence of 
the transport. IQtt et a l. (2012), using a neutrino leak- 
age scheme and stud ying the same 27Mq progenitor as 
IMuller et all (|2012a[ ). observed low- level SASI activity 
at early times, which was eventually suppressed in mod- 
els which exploded because of artificially enhanced neu- 
trino heating. Their neutrino treatment, however, is not 
on par with the mul ti-group ray-by-ray-plus transport of 
IMuller et all (|2012aD and therefore comparisons with the 
more sophisticated 2D models of IMuller et al.l (2012a) 
should be made only with great caution and reservation. 

In view of the poor exploration of conditions in col- 
lapsing stellar cores in 3D so far and considering the 
substantial approximations that have been made in 3D 
supernova simulations until now, the strong opinions ut- 
tered about the dominance of neutrino-driven convection 
(IBurrows et al.ll2012t [Murphy et al1l2012|; IDolence et al.l 
120131 ) and the categorical rejection of an impor tant role 
of the SASI in "realistic" 3D supernova models ([Burrowsl 
120131 ) are disturbing. Actually, it is not overly astonish- 
ing that the highly simplified setup investigated by these 
authors did not show signs of any pronounced SASI ac- 
tivity. The growth of SASI modes was clearly disfavored 
in their models by several aspects. Neglecting neutrino 
losses from the neutron star interior (above an optical 
depth of about unity) and using a relatively stiff nu- 
clear equation of state the authors prevented the neutron 
star from shrinking below ^60 km. Correspondingly, the 
shock radius remained rather large, in which case the 
postshock velocities were relatively small and the advec- 
tion timescale of the accretion flow through the gain layer 
was relatively long. This damped the development of the 
SASI, whose growth rate is roughly proportional to th e 
inverse of the advection timescale (cf. Schcck et al. 2008) . 
At the same time such conditions supported neutrino- 
driven convection, which preferably develops in the ac- 
cretion flow for ratios x °f the advection time to the 
local buoyancy times cale above a critical value of x ~ 3 
([Foglizzo et alJ[2006[ ). 

Although the theoretical understanding of the growth 
conditions of the SASI and neutrino-driven con- 
vection is mostly based on linear theory (e.g., 
Foglizzo et al.l [20071: lLaming||2007t lYamasaki fc Yamada 
20071 : lYamasaki fc Foglizzol 120081 iGuilet fe~ glizzo 



2012J), the predictions were found to be consistent 



with the behavior seen in 2D hydrodynamical simu- 
lations of accretion shocks in collapsing stellar cores 
(IScheck et al.l I2008D a nd in full-scale 2D supernova 
models (jMiiller et all l2012al) . Similarly, although 
it is a prior i not clear whethe r the SASI slosh- 
ing motions flBlondin et al.l 120031) and spir al modes 
([Blondin fc Mezzacappdl2007lrFe rnandcz 2010) observed 
in adiabatic accretion flows in 2D and 3D, respectively , 
or in shallow- water experiments ([Foglizzo et al.l [20121 . 
are conclusive for phenomena of relevance in the con- 
vectively unstable environment of the neutrino-heated 
layer, violent bipolar shock oscillations with SASI-typical 
characteristics were also identified to determine the evo- 
lution of the stalled supernova sho ck in some progen- 
itors and preexplos ion phases (e.g., IScheck et al.l 120081 : 
IMuller et al]|2012al) . While the phenomenon of the SASI 
as an advective-acoustic instability is not generically 
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linked to 2D, strong SASI activity has so far not been 
detected in full-scale 3D supernova simulations and it 
was speculated that its amplitude could be reduced by 
the kinetic energ y being shared betwe en three instead of 
two dimensions flwakami ct al. 2008), that the absence 
of a flow-constraining symmetry axis might disfavor co- 
herent mass m otions of low spherical harmonics modes 
( Burrows 2013), or that neutrino-driven buoyancy is gen- 
erally the fast est growing nonspher ical instability in su- 
pernova cores ([Burrows etaII[20T2] ). 

In this paper, we report about unambiguously iden- 
tified SASI activity in the first 3D simulation with de- 
tailed neutrino transport of the 27 Mp> progen itor model 
that was investigated by iMiiller et al.l (|2012a[) . For this 
3D simulation we employed the Prometheus- Vertex 
neutrino-hydrodynamics code with detailed ray-by-ray- 
plus neutrino transport and the sophisticated treatment 
of energy-dependent neutrino interactions al so applied in 
previous 2D simulations with this code (e.g. JBuras et al.l 
l2006bUMarek fc Jankd[2009l) a nd in the 2D GR models 
of IMiiller et al.l (|2012bi lal 120131) . Contrary to the claims 
discussed above, our model shows that despite the pres- 
ence of neutrino-driven convection, the SASI can grow 
no less vigorously in 3D (without any coordinate grid- 
imposed symmetry) than in 2D as long as small shock 
radii are maintained (in the 27 Mq star because of high 
mass accretion rates) and guarantee favorable growth 
conditions. We also observe the development of a clear 
spiral mode. SASI shock motions appear to be dimin- 
ished only when the accretion rate drops after the Si/SiO 
shell interface reaches the shock and the shock is able to 
expand to considerably larger radii. The variation of 
the relative strengths of neutrino-driven convection and 
SASI sloshing is consistent with experience and under- 
standing based on previous 2D simulations. In order to 
demonstrate this more directly we also performed para- 
metric 2D and 3D core-collapse and explosion simula- 
tions by using simp lified , gray neutrino transport as in 
Muller et all (|2012cft and iWongwathanarat et al.l (|2010bl 



20ll " to enhance the computational efficiency. Excising 



the high-density core and replacing it by a contracting 
inner boundary condition, at which neutrino luminosi- 
ties with chosen values could be imposed, allowed us to 
control the supernova core conditions. In these simula- 
tions we studied a 25 Mp, progenitor and made use of an 
axis-free Yin- Yang grid (jWongwathanarat et al.ll2010af ) 
instead of a polar grid. These simulations did not only 
show that the discussed SASI phenomenon can play a 
role in different progenitor stars and is not dependent on 
the choice of a particular computational grid, they also 
confirmed that SASI mass motions in 3D can be triggered 
by the same conditions that are known to be favorable 
in 2D. 

Our paper is structured as follows. In Sect. [2] we 
discuss the results of our fully self-consistent neutrino- 
hydrodynamics simulation of the 27 Mq progenitor, de- 
scribe briefly the numerical methods and modeling setup, 
analyze differences of the convective and SASI activity 
between 2D and 3D simulations, and present evidence 
for the temporary development of a spiral mode. In 
Sect. |3] we compare 2D and 3D results for three simu- 
lations of the 25 Mq model with our simplified modeling 
setup, whose main features we also summarize before we 
investigate the dependence of the development of SASI 



activity on the parametrically regulated conditions in the 
postshock and neutrino-heating layer. A discussion and 
conclusions will follow in Sect. [4] 

2. FULLY SELF-CONSISTENT 3D CORE-COLLAPSE 
SIMULATION OF A 27 Mq STAR 

2.1. Numerical methods and modeling setup 

The calculations for the 27 M© model in 2D and 
3D were performed with the elaborate neutrino- 
hydrodynamics code Prometheus- Vertex. This su- 
pernova simulation tool combines the hydr odynam- 
ics solver Prometheus (jFrvxell et al.l I1989T ) . which 
is a dimensionally-split impleme ntation of the piece- 
wise p arabolic method (PPM) of iColella fc Woodward! 
( 1984) , with the neutr ino transport module Vertex 
(Ram pp fe Jan ka 2002). Vertex solves the energy- 
dependent moment equations for the neutrino energy and 
momentum density (with full velocity dependence) for 
spherically symmetric transport problems defined to be 
associated with every angular bin of the polar grid ( "ra- 
dial rays") used for the multi-dimensional simulations. 
The moment equations are closed by a variable Edding- 
ton factor relation that is provided by the formal solu- 
tion of a model Boltzmann equation. An up-to-date set 
of n eutrino interaction r ates is included in Vertex (see, 
e.g.. IMiiller et al. 2012b). In the multi-dimensional case, 
our ray-by-ray-plus approach (jBuras et al.l l2006bl ) in- 
cludes non-radial neutrino advection and pressure terms 
in addition to the radial transport solves. In the sim- 
ulations presented here, we assume monopolar gravity, 
but include general relativistic cor rections by means of 
an effective gravitational potential ()Marek et al.l l2006) . 

We simulate the e volution of the 27M© progenitor of 
iWooslev et al.l (120021) . w hich was previously inv estigated 
bv IMiiller et all (|2012al ) and lOtt et all (|2012f ). both in 
2D and i n 3D, using the hig h -dens ity equation of state 
(EoS) of Latt imer fc Swestvl (|1991f ) with a nuclear in- 
comprcssibility of K = 220 MeV. The models are com- 
puted on a spherical polar coordinate grid with an initial 
resolution of n r x ng x n v = 400 x 88 x 176 (3D) and 
n r xng — 400 x 88 (2D) zones. Later, refinements of the 
radial grid ensure adequate resolution in the PNS surface 
region. The innermost 10 km are computed in spherical 
symmetry in both cases to avoid excessive time-step lim- 
itations. Seed perturbations for aspherical instabilities 
are imposed by hand 10 ms after bounce by introducing 
random perturbations of 0.1% in density on the whole 
computational grid. 

2.2. SASI activity: 2D versus 3D 

While both the SASI and convection can lead to large- 
scale shock deformations, the SASI is distinguished by 
a characteristic oscillatory growth and in its nonlinear 
stage by the quasi-periodic, oscillatory nature of the 
shock motions. In 2D, the artificial symmetry constraint 
and the excitation of large-scale modes by the inverse 
turbulent cascade could still produce a quasi-periodic 
sloshing motion even in convectively-dominated models 
([Burrows et alll2012l ) , and a more refined analysi s is nec- 
essary to identify the SASI (|Muller et alJl2012aft . How- 
ever, in 3D the distinction is much clearer, since large- 
scale shock deformations caused by buoyancy-driven con- 
vection initially evolve randomly without any identifiable 
periodicity and then grow monotonically once they reach 
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Figure 1. Snapshots of phases with convective and SASI activity in the evolution of the 27 Mq model at 154 ms, 223 ms, 240 ms (upper 
panels, from left to right), 245ms, 249ms, and 278ms (lower panels, from left to right). The volume rendering visualizes surfaces of 
constant entropy: The outer, bluish, semi-transparent surface is the supernova shock, the red surfaces are entropy structures in the 
postshock region. The upper left panel displays mushroom-like plumes of expanding, high-entropy matter that are typical of neutrino- 
driven buoyancy. The upper middle and right plots and the lower left and middle panels show distinctly different entropy structures 
of dipolar (and quadrupolar) asymmetry, which engulf the still visible buoyant plumes with their higher-order spherical harmonics mode 
pattern. The entropy asymmetries of I = 1, 2 character are caused by global shock sloshing motions, which create hemispheric high-entropy 
shells in phases of shock expansion. At 223 ms and 240 ms the shock has pushed towards the lower right corner of the panels whereas 
at 245 ms and 249 ms it is in a phase of violent expansion motion towards the upper left corner of the plots. All stages exhibit a strong 
deformation of the shock. At 278 ms the vivid SASI phase is over, the shock is more spherical again, and the postshock entropy structures 
correspond to neutrino-driven plumes. 



a certain threshold amplitude ([Burrows eTaII[20Tl . Pe- 
riodic SASI oscillations and large-scale shock deforma- 
tions caused by convection can therefore hardly be mis- 
taken for each other in 3D. 

Images of the entropy distribution in the postshock 
layer (Figure [T]) and in particular the corresponding 
movie of our 3D simulation of the 27 Mq progenitor 
indeed provide a clear hint that both distinctly differ- 
ent instabilities are at work in the shocked accretion 
flow around the nascent neutron star. The instabili- 
ties develop nearly at the same time and are present si- 
multaneously for an extended period of the simulated 
postbounce evolution. The first small mushroom-like 
Ray leigh- Taylor fingers of neutrino-driven convection be- 
come visible around 80-100 ms after bounce to subse- 
quently grow stronger and larger in angular size over a 
timescale of some ten milliseconds. At about 125 ms p.b. 
the rising plumes begin to cause shock deformation and 
a modest amount of global asphericity of the accretion 
layer. Until ~155 ms the activity in the postshock layer 
is clearly dominated by neutrino-driven buoyancy (Fig- 
ure [H upper left panel), but at t > 155 ms, during a phase 



of accelerated shock recession, coherent entropy struc- 
tures show up first. The corresponding low-mode spher- 
ical harmonics pattern clearly differs from the buoyant 
mushrooms on smaller angular scales. This phenomenon 
is associated with shock oscillations, which quickly am- 
plify to bipolar shock sloshing motions and create charac- 
teristic, hemispheric high-entropy shells during phases of 
fast shock expansion. These half-shells of shock-heated 
matter engulf the buoyant bubbles of neutrino-driven 
convection in deeper regions (Figure [TJ upper middle and 
right and lower left and middle panels). While the slosh- 
ing axis initially wanders, it becomes more stable as the 
SASI sloshing of the shock further grows in amplitude 
and violence between ^195 ms and ^240 ms. As a con- 
sequence, an expansion of the average shock radius is 
driven even before the Si/SiO composition-shell interface 
arrives at the shock and the mass accretion rate starts to 
drop steeply at t ~ 220ms (Figure^ upper left panel). 
The decrease of the accretion rate supports the shock 
expansion, in course of which the bipolar, quasi-periodic 
shock pulsations gain even more power. At t ~ 225 ms 
a spiral mode seems to set in for several revolutions be- 
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Figure 2. Comparison of SASI activity and shock evolution in the 2D and 3D simulations of the 27 Mq model. Panel 
radius ((r^) = ao) an d mass accretion rate of the collapsing stellar core at 400km; panel b: components of the SASI 
vector; panel c: rms shock deformation <r; panel d: kinetic energy of non-radial mass motions in the gain layer. 
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Figure 3. Power spectra of v v of the 27 Mq 3D model sam- 
pled between r = 63 km and 80 km at post-bounce times of 222 ms 
(black) and 259 ms (red). At 222 ms, the strong SASI produces 
a distinctive peak at I = 1, which is absent during the later, 
convection-dominated phase. 



fore the average shock radius reaches its maximum ex- 
tension at ~ 250 ms and the SASI sloshing dies off at 
t > 260 ms. The presence of large-amplitude spiral mo- 
tions is reflected by considerable variations of the mean 
shock radius between 230 ms and 260 ms. These disap- 
pear when the SASI activity ceases at t ~ 260 ms (Fig- 



ure [H upper left panel). Later on, until the end of our 
3D simulation, aspherical mass motions in the postshock 
layer are dominated again by the buoyant plumes typi- 
cal of neutrino-driven convection (Figure [TJ lower right 
panel) . 

This verbal description of the dynamical evolution 
of the postshock accretion layer is supported by a de- 
tailed analysis based on several timc-dependent parame- 
ters that quantify the characteristic features of SASI ac- 
tivity. To this end we perform a time-dependent decom- 
position of the angle-dependent shock position r s h(#, </?) 
into spherical harmonics Y™: 



(-1)' 



»l 



y/4ir(2£+ 1) 



r sh (9,v)Y e m (e,p)dn. (1) 



Here the Y™ are real sphe rical harmonics with the 
same normalizat ion as used bv lBurrows et all ()2012[ ) and 
lOtt et all (|2012D . With this choice of basis functions, the 
coefficients with 1 = 1 give the angle-averaged Cartesian 
coordinates of the shock surface, 



(Vsh) 



Hi' 



(Zsh) 



(2) 

and a[j is just the average shock-radius (r s h). 

The time evolution of the coefficients a x , a y (3D), and 
a z (3D and 2D) is shown in panel b of Figure [2] Both 
in 2D and in 3D, the shock surface clearly oscillates in 
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Figure 4. Evolution of the stability parameter \ for the gain layer 
of the 3D simulation of the 27 Mq progenitor. During most of the 
time \ < 3. This suggests conditions in the postshock accretion 
flow which disfavor the growth of neutrino-driven convection rela- 
tive to the d evelopment of the S ASI in analogy to the 2D model 
discussed bv lMuller e"t""aH H2012al) . 



a quasi-periodic manner until ~260ms after bounce, i.e. 
until shortly after the Si/SiO shell interface has reached 
the shock and the accretion rate has dropped consider- 
ably between ~220ms and 240 ms p.b. The lower ac- 
cretion rate results in a pronounced expansion of the 
average shock radius (panel a of Figure [2]) , which ini- 
tially is stronger in 3D. However, the 2D model maintains 
large (albeit less regular) shock oscillations, with the av- 
erage shock radius eventually overtaking the 3D model 
at ~ 300 ms when an explosion de velops (i.e. somewh at 
later than in the GR simulation of lMuller et al.ll2012aTFl . 
By contrast, the shock continues to recede in the 3D run. 
The more optimistic evolution of the 2D model compared 
to the failing 3 D model at late sta ges is consistent with 
the findings of lHanke et al.l (|2012l ). and could be due to 
the action of the inverse turbulent energy cascade, which 
continues to feed energy into large-scale modes in 2D. 

However, the evolution of the two simulations prior to 
the infall of the Si/SiO interface is remarkable: While the 
amplitude of the t = 1 mode is initially larger in 2D, the 
individual components a x , a v , and a z of the 1=1 ampli- 
tude vector in 3D become comparable to a z in 2D around 
200 ms, and a x even reaches considerably bigger values. 
During this phase, the SASI is undoubtedly stronger in 
3D than in 2D. Further confirmation of this assessment is 
provided by the root-mean-square deviation cr(r s h) of the 
shock radius from its average value (panel c of Figure [2]) : 



(47T) 



(r s h(0,^)-(r s h» 2 dfi 



(3) 



For reasonably small amplitudes, a is also a measure for 

1 It is not clear whether this difference or how much of this dif- 
ference is caused by GR effects, because the models in the present 
paper were simulated with a slightly different treatment of the low- 
density equation of state, which led to a significant delay (~35 ms) 
of the infall of the silicon layer and a correspondingly later arrival 
of the Si/SiO shell interface at the shock. 
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Figure 5. Snapshot of the entropy (color coded according to color 
bar at the upper left corner, in units of Boltzmann's constant fcj, 
per nucleon) in the plane through the origin normal to the vec- 
tor n = (—0.35,0.93,0.12) at a post-bounce time of 152 ms in the 
27 Mq 3D model. The high-entropy plumes with high-order spheri- 
cal harmonics pattern suggest buoyancy-driven convective overturn 
of neutrino-heated matter. 



the total power of SASI amplitudes with different i: 



\ 



EE 

1=1 m=-l 



(4) 



The same picture emerges when we consider the kinetic 
energies Eg and E v associated with motions in the 9- 
and (^-directions in the gain region, 

E e = l[ PWfldV, E v = \ 

J Vg-ain 



pv% dV. 



(5) 



As shown in panel d of Figure [3J the total energy con- 
tained in non-radial motions is also larger in 3D during 
the relevant phase around ^230 ms. In the period of con- 
tinuous increase of the SASI amplitude in the 3D model 
between t ~ 155 ms and ~ 240 ms, the kinetic energy 
grows and a exhibits quasi-periodic modulations signal- 
ing the shock sloshing motions. Interestingly, during the 
phase of strongest SASI activity we find rough equipar- 
tition between the kinetic energies of non-radial motions 
Eg + E v , and the energy E r contained in fluctuating ra- 
dial velocities, 

1 



E r 



p(v r -(v r }) dV, 



(6) 



where (v r ) is the angle-averaged radial velocity. This 
cquipartition is apparently no t a unique feature of 
buoyancy-driven turbulence (cf. iMurphv et al.ll2012l) . at 
least not as far as these volume-integrated quantities are 
concerned. 

A clear difference between SASI dominated and con- 
vection dominated phases of the 27 Mq 3D model can 
be observed in the power spectrum of the azimuthal ve- 
locity v v as a function of multipole order I. Figure [3] 
shows the spectra during a SASI active phase (222 ms 
p.b.) compared to the later time (259 ms) when the SASI 
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motions cease and convective plumes with their higher- 
order multipole pattern determine the asphericities in the 
postshock region ag ain. The pow e r spec tra are evaluated 
with equation (8) of Hankc et al.l (|2012l ) for v v (weighted 
with the square root of the density) integrated over a 
radial region between 63 km and 80 km. The presence 
of the SASI low- mode deformation at 222 ms leads to a 
prominent peak of the power spectrum at low multipole 
orders, which is absent in the later spectrum. 

The shock deformation and SASI amplitude moder- 
ately increase as the shock retreats between ^100 ms and 
200 ms p.b. (The same is true for the specific nonra- 
dial kinetic energy in the gain layer, although the kinetic 
energy decreases temporarily because of the decreasing 
mass in the gain region.) It is noteworthy that the few 
tens of milliseconds of stronger SASI activity in 3D coin- 
cide with a phase of more rapid shock expansion than 
in 2D. Conceivably, the more energetic SASI motions 
provide a stronger push against the pre-shock ram pres- 
sure. A close examination of the (r s h) m Figure [2] shows 
that some (oscillatorily modulated) SASI-aided shock ex- 
pansion seems to set in around 190 ms, i.e. already be- 
fore the rapid drop of the preshock mass-accretion rate 
that begins at ~220ms. Therefore one might speculate 
that with slightly more time available for the growth of 
the SASI, the extra support by nonradial SASI motions 
might have driven the 3D model over the threshold for a 
neutrino-powered runaway expansion of the shock after 
the infall of the Si/SiO interface. 

Quite remarkably, the SASI is not only able to reach 
larger amplitudes in 3D than in 2D as long as its growth 
conditions remain favorable, but it is even found to 
develop despite some earlier convective activity. Fig- 
ure |4] displays the critical parameter y for the growth of 
convection as evaluated from spherically averaged stel- 
lar quantities in the gain layer of our 3D simulation 
of the 27M( 7) proge nitor according to equation (3) in 
IMiiller et alf p012ah . 

(7, 

J(r e ) I ( V r) 

where wbv is the Brunt- Vaisala frequency. The integra- 
tion is performed between the average gain radius (r g ) 
and the average shock radius (r s h). Note that only re- 
gions contribute to the integral where Wg V < indicates 
local instability. The parameter y roughly measures the 
ratio between the advection timescale of the flow through 
the gain layer and the growth timescale of convection. 
Since perturbations are advected out of the gain layer 
with the accretion flow in a finite time, convection can de- 
velop only when perturbations are amplified sufficiently 
strongly within this time. F or the linear regime ( i.e., for 
small initial perturbations) Foglizzo et al.1 (p006| ) found 
the threshold condition of y > 3 for convective activity 
to develop in the accretion flow of the gain layer. This 
result of mathematic al analysis is suppor t ed by numer- 
ical st u dies in 2D bv iBuras et al.l (l2006af) ; IS check et al.l 
(2008ft: lFer"nandez fc Thompson! (|2009bl laj). 

Despite y < 3 (FigureSJ), however, convection develops 
around 80 ms after bounce in our 3D simulation of the 
27 Mq model. This happens because convective activity 
is not only seeded by the artificially imposed, random 
density perturbations of 0.1% amplitude (cf. Sect. 12. ip 



but also by numerical perturbations along the axis of 
the computational polar grid in one hemisphere, which 
we are not able to damp perfectly. Although still small, 
these numerical effects are sufficiently large to trigger 
the rise of a buoyant plume against the advection flow, 
which instigates further perturbations that exceed the 
linear regime. In this case convection can be initiated 
altho ugh y < 3 signals st ability according to linear anal- 
ysis (|Foglizzo et al.l [2006ft . We note in passing that the 
level of seed perturbations was smaller and well compat- 
ible wit h linear theory in the 2D models of IMiiller et al.l 
(|2012a|) . and we also emphasize that the axial perturba- 
tions have a noticeable effect only in the early growth 
phase of convection whereas no significant axial artifacts 
or alarming flow perturbations near the polar axis can be 
observed during the later phases of fully developed non- 
radial flow activity in the postshock flow (see Figure [lj . 
Once strong SASI and convective mass motions have de- 
veloped in the flow, the supernova core contains a noise 
level so high that a subsequent decline of y below the 
critical threshold of ~3 does not imply that conv ective 
activity is unable to continue (jMuller et al.ll2012aft . 

Figure [5] shows a snapshot of the entropy in a 2D slice 
at 152 ms. Here, the post-shock flow is still dominated 
by multiple, intermediate-scale plumes as familiar from 
buoyancy-driven convection. SASI shock sloshing be- 
comes strong and temporarily dominant only afterward. 
The argument that any convective activity arising from 
sufficiently large seed perturbations will quench the SASI 
thus seems to be invalid. The competition between con- 
vection and the SASI is obviously more subtle than a 
superficial reading of recent papers dBurrows et al.ll2~012l ; 
IMurphv et all 12012: IMiiller et alJl2012ah IQtt et al.ll20T2ft 
might suggest. Some relevant aspects of the competing 
growth conditions and interaction of the two instabilities 
we re discussed on the basis of 2D supernova simulations 
bv lScheck etafl (12008ft . A very similar behavior can be 
diagnosed in the 3D case and will be addressed on the 
basis of parametric 3D studies in Sect. |3l where we will 
make use of an axis-free Yin- Yang grid in order to avoid 
the perturbative influence of the polar coordinate axis 
during the linear growth phase of the seed perturbations. 

2.3. Detecting the spiral mode of the SASI 

While the SASI is limited to a sloshing motion along 
the symmetry axis in 2 D, there is the possibi l ity of 
a spiral mode in 3D (iBlondin fc Mczzaca ppal 120071 : 
llwakami et all 2009; Fcrnandeil2010), which could pro- 
vide a means for angular momentum separation between 
the PNS and the ejecta, and might also have different 
saturation properties in the non-linear phase. 

Detecting the spiral mode is not straightforward, as 
the time evolution of the coefficients a™ needs to be 
taken into account. Merely computing the coefficients 
with m ^ is not sufficient. In principle, spiral and 
sloshing modes can be disenta ngled by a Fourier analysis 
of a™ (t) (jlwakami et al.ll2008D if they remain stable over 
several oscillation periods. We use a somewhat different 
approach to visualize the character of the I = 1 mode in 
our 3D model here. The coefficients a™ can be combined 
into a vector, 

ai = (a^af 1 , a?) = (a^a^a*), (8) 

which is a rough measure of the angle-averaged displace- 
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Figure 6. Evolution of the i = 1 amplitude vector ai for the 27 Mq simulation from two different viewing angles. Different colors are 
used for the phase up to 177 ms (yellow), the phase of strong SASI activity (blue, up to 265 ms), and the late, SASI-quiet phase (red). 



merit of the shock center from the origin and also indi- 
cates the direction and amplitude of the shock deforma- 
tion (neglecting modes with higher £). We visualize the 
time evolution of this amplitude vector in 3D space in 
Figure [5] with different colors of the trajectory indicating 
different phases of the simulation. 

Once the SASI starts to grow vigorously (blue curve), 
ai initially moves along a narrow elliptical path with 
growing semi-major axis, indicating a predominant slosh- 
ing mode. Towards the phase of strongest SASI activity, 
the trajectory becomes more circular, signaling the tran- 
sition to a spiral mode. The plane of the spiral remains 
relatively stable until the maximum amplitude is reached 
and the SASI dies down again. It is roughly perpendic- 
ular to the vector n = (—0.35, 0.93, 0.12), i.e. there is no 
alignment with the axis of the spherical polar grid. Fig- 
ure [S] also further illustrates the different behavior during 
the "SASI-dominated" phase compared to the earlier and 
later "convectively-dominated" phases, during which ai 
evolves in a more random fashion. 

Slicing the model along the plane in which ai predom- 
inantly moves allows us to visualize the distinctive spi- 
ral mode pattern, as shown in Figure [7] The snapshots 
of the rotational velocity around the "axis" of the spi- 
ral mode (in the plane through the origin perpendicular 
to that axis) reveal two counter-rotating regions. While 
these regions are initially of comparable size, the flow 
in the counter-clockwise direction eventually dominates, 
and the rotation of the mode pattern with a continuously 
shifting triple point can clearly be seen. 

During the short phase of strong SASI activity, angu- 
lar momentum separation by the spiral mode proceeds 
very efficiently, transferring a total angular momentum 
of ~ 5 x 10 46 erg s into the gain region (Figured]). As 
expected, the direction of the angular momentum vector 
is extremely similar to the normal vector of the spiral 
plane. However, all the angular momentum is eventually 
advected out of the gain region after the SASI has died 
down. 



3. PARAMETRIC 2D AND 3D POST-BOUNCE 
SIMULATIONS OF A 25 M© STAR 

In order to probe the SASI growth conditions and 
the differences between 2D and 3D accretion flows in 
more detail, we performed a set of s imulations for 
the 25 Mp~, solar- metallicity progenitor of iWooslev et al.l 
(2002). The choice of this stellar model is motivated by 
the fact that it turned out to be particularly hard to 
explode in self-consistent 2D simulations with sophisti- 
cated neutrino treatment, where it exhibits violent shock 
sloshing motions with amazingly stab le periodicity over a 
long period of po stbounce evolution (jMiiller et al.ll2013l : 
Uanka et al.ll2012l ). It thus offers promising perspectives 
to explore the relevance of the SASI in 3D for a second 
progenitor in addition to the 27 Mq case investigated in 
Sect. H 

3.1. Numerical methods and modeling setup 

In order to efficiently run several 3D simula- 
tions, a simplified setu p is employed that was exten- 
sively applied befor e by iWongwathanarat et al.1 (l2010bl 
20121): iMuller et alJ (l2012cD: iScheck et all (120061 120081k 
Arcones et al.l (|2007| ): lArcones fc Jankal (|2011l ). Instead 
of the elaborate neutrino transport solver used for 
the 27 Mq run in Sect. [21 now an approximate, gray 
(ray- by-ray ) neutr ino transport is used (for details, see 
IScheck et al.ll200"6T) . Moreover, the high-density interior 
of the proto-neutron star (at densities above neutrino 
optical depths between some 10 and several 100, de- 
pending on the phase of the evolution) is excised and re- 
placed by an inner grid boundary (also in contrast to the 
27 Mq run, where the whole neutron star was included 
in the simulation). A contraction of the boundary radius 
is prescribed to mimic the shrinking ne utron-star core 
(jScheck et al.l 12006b lArcones et all 120071 ). and neutrino 
luminosities of chosen magnitude and time dependence 
are imposed there to parametrize the neutrino losses from 
the interior core volume. The transport approximation 
allows us to account for radial flux variations due to ac- 
cretion luminosity and neutrino energy deposition in the 
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Figure 7. Snapshots of the rotational velocity around the origin in the plane perpendicular to n = (—0.35,0.93,0.12) at post-bounce 
times of 223 ms, 227 ms, 231ms, and 235 ms during the 3D simulation of the 27 Mq progenitor. Red and yellow (positive velocity values) 
correspond to counterclockwise rotation. 
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Figure 8. Time evolution of the angular momentum components 
L x , L y , and L z , and of the absolute magnitude \L\ of the total 
angular momentum contained in the gain region of the 27 Mq sim- 
ulation in 3D. 



gain layer. Although this modeling approach is not elab- 
orate enough to exactly reproduce the results obtained 



with more sophisticated and fully consistent treatments, 
it nevertheless includes a reasonably good representation 
of all aspects that determine the flow dynamics in collaps- 
ing stellar cores. The parametrized approach has the big 
advantage that different components can be controlled to 
a large extent independently so that studies with varied 
conditions can be performed in a systematic manner. 

The computational efficiency is further en- 
hanced by the use of an axis-free Yin- Yang grid 
(jWongwathanarat et al.l l2010af l . which avoids the time- 
step constraints and numerical artifacts close to the 
axis of the polar grid. The 3D runs were conducted 
with a grid zoning of 600(r) x 47(0) x 137 (tp) x 2, 
the corresponding 2D models with 600(r) x 90(6) 
zones, which ensures an angular resolution of 2° in 
both 2D and 3D. The inner grid boundary is placed 
at an enclosed mass of approximately 1 . f Mq . When 
the models are started from ID conditions at 15 ms 
after bounce, the corresponding initial radius of the 
inner grid boundary, R\ h , is located at 53 km, while 
the radius of the outer grid boundary, i? Q b, is chosen 
to be at 10000 km in the beginning. The radial grid 
resolution, Ar, is kept constant at 0.265 km up to the 
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radius of initially 100 km, beyond which the zone width 
is increased logarithmically. 

The shrinking core of the nascent neutron star is de- 
scribed by the contraction of the inner grid b oundary 
according to equation (1) of Scheck ct al. (2006]) with an 
exponential timescale of tib = 1 s and a final radius of 
R\ h . Along with the retraction of the inner grid bound- 
ary the whole radial grid is moved with the same velocity. 
During the computed model evolution the neutrino lu- 
minosities imposed at the inner grid boundary are kept 
constant. The mean neutrino energies of the neutrino 
fluxes entering the grid are prescribed as functions of the 
temperature in the first radial cell of the computational 
grid as described bv lUgliano et al.l (|2012D . 

In order to trigger the growth of nonradial hydrody- 
namic instabilities, seed perturbations with randomly 
varied amplitude are imposed as zone-to-zone variations 
of the radial velocity in the whole computational domain 
at the start of the 2D or 3D simulations. 

We present results of three types of models, which dif- 
fer in the size of the seed perturbations, either of 0.1% or 
3% amplitude ("small" and "large" initial perturbations, 
respectively), and in the contraction behavior adopted 
for the neutron-star core. The latter is controlled by the 
assumed asymptotic radius of the inner grid boundary 
at t — > oo, namely R ib = 15 km for "slow contraction" 
and R l ih = 10km for "fast contraction". In the case of 
slow contraction the imposed boundary luminosities of v e 
and v e are 1.49 x 10 51 ergs -1 and 0.89 x 10 51 ergs -1 , re- 
spectively; in the case of fast contraction we use L v ^ jb = 
5.96x10 s1 ergs -1 and = 3.58x 10 51 ergs -1 . (Muon 

and tau neutrinos and antincutrinos are included as well 
but do not play an important role for the hydrodynamics 
of the accretion flow in the supernova core.) The three 
computed models are characterized by the following pa- 
rameter sets: 

S25-1: slow boundary contraction, 3% perturbation am- 
plitude; 

S25-2: fast boundary contraction, 0.1% perturbation 
amplitude; 

S25-3: fast boundary contraction, 3% perturbation am- 
plitude. 

Model S25-1 was evolved until 515 ms after bounce, 
while S25-2 and S25-3 were run until 365 ms after core 
bounce. 

3.2. SASI activity in dependence of supernova-core 
conditions 

Figure [9] provides an overview of the dynamical evolu- 
tion of the three models S25-1, S25-2, and S25-3 in 2D 
compared to 3D by showing the north polar and south 
polar entropy profiles of the 2D simulations in the top 
panels, shock radii (average, maximum, minimum) in the 
second line from the top, the critical x parameter (eval- 
uated according to Equation [7]) in the panels of the third 
line, and a spherical harmonics mode decomposition of 
the deformed shock surface in the fourth and fifth lines 
for 2D and 3D cases, respectively. 

Here we follow the defo r matio n mode analysis used 
by IWongwathanara t et al.l (|2012f) and write the time- 
dependent decomposition of the 3D shock surface, 



r s h(#, if), in spherical harmonics as 



(9) 



where &™ are the expansion coefficients and the spherical 
harmonics Yf 1 are functions of the associated Legendre 
polynomials PJ™, 



Y e m (o, ip) = Kppp (cos oy™* . 



with 
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Multiplying Equation © by the complex conjugate of 
the spherical harmonic, Y" 1 *, and integrating over the 
solid angle, the expansion coefficient is found to be 

bf=l dpi d6 sin 9 r sh (6 , ip)Y e m * (9 ,cp). (12) 
Jo Jo 

For our mode analysis we actually use the pseudo-power 
coefficients c\ = |6q| 2 and 
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for £ > 0, respectively, whose values do not depend on 
the particular choice of the orientation of the coordinate 
grid in 3D. The cj are related to the coefficients a™ of the 
real spherical harmonics decomposition in Equation (QJ 

by 
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(14) 



In Fi gure [^] logarithmic values of c 2 normalized by Cq are 
plotted. 

The north polar and south polar entropy profiles of the 
2D simulations in the upper panels of Figure [9] show that 
the slow contraction case, S25-1, does not explode within 
500 ms after bounce, whereas the fast contraction mod- 
els S25-2 and S25-3 develop explosions around 300 ms 
post bounce because of the stronger neutrino heating by 
higher boundary and, in particular, accretion luminosi- 
ties. The dynamics of the shock and postshock flow until 
the explosions set in reveals interesting differences in all 
three cases. 

The north and south polar shock trajectories as well as 
the maximum and minimum shock radii (given by dot- 
ted lines in the panels on the second line of Figure [9]) 
in all three 2D cases exhibit pronounced sloshing mo- 
tions, which suggest considerable activity on the scale 
of low-order spherical harmonics modes. The most reg- 
ular pattern of alternating north-polar and south-polar 
shock expansions and contractions, however, is visible 
in model S25-2, whose periodic shock excursions remind 
one strongly of the clean SASI sloshing mode found for 
the 27 M Q progenitor by IMiiller et all (|2012aD . Before 
large-amplitude shock motions are reached, the slosh- 
ing amplitude grows slowly in an otherwise seemingly 
quiet environment. A similar behavior was observed by 
iScheck et all (j2008) in their 2D models WOO and WOOF, 
whose parameter settings (chosen core contraction rate 
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Figure 9. Comparison of the time evolution of 2D and 3D results of our set of parametrized postbounce simulations S25-1, S25-2, and 
S25-3 of the 25 Mq progenitor. The upper panels show north polar and south polar entropy profiles for the 2D runs with thin white lines 
indicating the shock trajectories, which separate low-entropy preshock matter (black and blue) from the high-entropy postshock region 
(red, orange, yellow). The second panels from the the top display maximum, average, and minimum shock radii for 2D (dotted) and 3D 
models, the panels in the third line present the \ parameter for convective instability (Equation of 2D (black lin es) a nd 3D (red lines) 
cases, and the plots in the fourth and fifth lines provide the normalized pseudo-power coefficients c?/cq (Equation 1 13 6 on a logarithmic 
scale for the 2D and 3D simulations, respectively. Phases with values of \ < 3 correlate with preferred growth of low-mode (I = 1, 2) SASI 
activity. 
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Figure 10. Structure of the 2D (left) and 3D (right) models of our set of parametrized 25 Mq simulations (top: S25-1; middle: S25-2; 
bottom: S25-3) at representative postbounce times (as given in the upper right corner of each panel). The central object in all panels 
provides a three-dimensional visualization of the supernova shock (white, transparent surface), which engulfs a surface of constant entropy 
of 17fcb P er nucleon. The images on the cube walls show entropy distributions (the color coding as given by the color bars in the lower left 
corners with black, blue, and green signaling low values) in the x-y, x-z, and y-z planes (orientation according to triad in the lower right 
corner of the upper left panel). The white circle indicates the inner grid boundary and a scale stick in each panel gives a measure of the 
size of the displayed volume. While the 3D model of S25-1 exhibits signatures of neutrino-driven convection and of low-order multipolc 
asymmetry due to SASI shock sloshing, both the 2D and 3D cases of S25-2 show a clear dominance of SASI activity, whereas model S25-3 
possesses the clearest pattern of neutrino-driven buoyancy. 
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in combination with very low core-boundary luminosi- 
ties) were initially disfavorable for the onset of neutrino- 
driven convection and thus allowed one to follow the 
SASI growth from the linear to the nonlinear regime in 
a very clear way. 

In contrast, the 2D models of S25-1 and S25-3 show 
shorter-timescale variability of the shock trajectories su- 
perimposed on a longer-timescale quasiperiodic north- 
south sloshing of the shock. This suggests the presence 
of asymmetries of higher-order spherical harmonics type 
as they are typical of neutrino-driven buoyant plumes. 

These descriptions are supported by a more quantita- 
tive analysis of the growth conditions of convective insta- 
bility in the postshock accretion layer as expressed by the 
critical x parameter, which measures the ratio of aver- 
age advection timescale of the accretion flow through the 
gain region to the convective growth timescale. Inspect- 
ing the panels in the third line of Figure [SJ where the 
black curves display the 2D results, reveals that only in 
model S25-2 \ drops below the critical threshold value of 
3 for most of the postbounce evolution until strong shock 
expansion and explosion sets in. The low values of x are 
a consequence of the strong shock contraction in response 
to the rapidly shrinking core of the neutron star assumed 
in this mo del. As expected from the linear perturbation 
analysis of Foglizzo ct al. (2006), which is applicable to 
model S25-2 because of the small amplitude of the initial 
seed perturbations, convective overturn cannot develop 
and the model exhibits the observed clean SASI growth 
pattern. Correspondingly, the spherical harmonics de- 
composition of the shock deformation shows the charac- 
teristic oscillatory growth of dominant £ = 1 (dipole) and 
1 = 1 (quadrupole) modes to nonlinear magnitudes over 
a timescale of roughly 100 ms (fourth panel from top in 
the middle column of Figure [9]) . 

In contrast, the 2D models of S25-1 and S25-3 show a 
growth of the shock asymmetries that differs from that 
of S25-2. A long-time oscillatory rise of the amplitude is 
not visible in both cases, but relatively large shock defor- 
mation develops much earlier and the time dependence of 
the spherical harmonics amplitudes of order i = 1, 2, 3 is 
much more irregular than in S25-2. Although still dipolar 
and quadrupolar shock deformation modes dominate, the 
short-timescale variability suggests bubbles and plumes 
on the smaller angular scales typically connected to the 
occurrence of strong buoyancy activity in the neutrino- 
heating layer. This is in agreement not only with x va l~ 
ues around the critical number of 3 for S25-1 and values 
well above 3 for S25-3, but it is also compatible with 
the visual impression of the entropy structures between 
shock and gain radius displayed in the left panels of Fig- 
ure 1101 The pronounced deformation along the symme- 
try axis, which is responsible for the dominance of the 
dipolar mode in many evolution phases of all three mod- 
els, can be nicely seen in the shape of the shock and 
of the surface of constant entropy displayed by the cen- 
tral, three-dimensional representation of the 2D model 
structure in each of the left panels. It should be noted, 
however, that different from the situation of S25-1 and 
S25-3, where the entropy surface shows short-wavelength 
variations reflecting the presence of buoyant plumes, this 
surface is much smoother in model S25-2, where the vio- 
lent SASI sloshing at 280 ms post bounce is accompanied 
by much weaker convection. 



The rapid onset of buoyancy in models S25-1 and S25-3 
is fostered by the larger initial seed perturbations. Since 
the difference between S25-2 and S25-3 is only the size 
of the amplitude of the initial random seeds (0.1% for 
S25-2 and 3% for S25-3), the different evolution of the 
flow asymmetries and the different values of x m these 
two models clearly demonstrate the important influence 
of the magnitude of preexisting deviations from spherical 
symmetry in the collapsing stellar core. 

A comparison of 2D and 3D results for the same mod- 
eling setup is extremely interesting now. Inspecting the 
right panels of Figure [10] it becomes obvious that the 
SASI activity appears in a much clearer way in the en- 
tropy surfaces of the 3D models. The shock sloshing 
creates smooth, coherent hemispheric shells of constant 
entropy like those we have seen in Figure [1] and al- 
ready discussed in Sect. 12.21 The 3D case of S25-2 ex- 
hibits such bipolar hemispheric entropy lobes, which are 
nearly axisymmetric around a direction tilted relative 
to the three coordinate axes (indicated by the triad), 
but which possess hardly any substructure on shorter 
wavelength scales. This indicates that the entropy vari- 
ations are created by pure SASI sloshing movements of 
the shock. The three x-y, x-z, and y-z cuts displayed 
on the cube walls yield an impressive confirmation of 
nearly perfect entropy half-shells without any indication 
of secondary, parasitic instabilities of Rayleigh- Taylor or 
Kelvin-Helmholtz type. In contrast, model S25-3 shows 
the familiar cauliflower pattern of expanding and frag- 
menting Rayleigh- Taylor fingers and Kelvin-Helmholtz 
swirls typical of neutrino-driven buoyancy. Model S25-1 
is again different and characterized by features of both 
types in superposition: The snapshot of the upper right 
panel in Figure [10] displays a convective bubble pattern 
that is engulfed by a smooth half-shell closer behind the 
shock in one hemisphere. 

These differences are reflected by the spherical har- 
monics amplitudes for the I = 1, 2, 3 modes of the shock 
deformation in the 3D models presented in the bottom 
panels of Figure M While in models S25-1 and S25-2 
the I = 1 deformation around the times of the snap- 
shots in Figure [T0l (330 ms and 280 ms, respectively) re- 
mains clearly dominant for longer evolution phases, and 
the time variation of the dipole amplitude reveals very 
stable periodicity for 100 ms and more, all three lowest- 
order spherical harmonics have similar strength in model 
S25-3 around 280 ms after bounce. In S25-3 the dipolar 
asymmetry component has the largest amplitude only 
for some short episodes of the postbounce evolution, al- 
though quasi-periodic 1=1 amplitude variations can 
be observed also in this case. The higher-order multi- 
poles in all models show greater variation frequencies and 
less good temporal regularity and seem to be influenced 
by the stochasticity associated with convective mass mo- 
tions. 

The x parameter (panels in the third line of Figure [§]) 
aids a better understanding of the differences between 
the 3D models and of their differences relative to the 
corresponding 2D results. In the 2D and 3D cases of 
model S25-2 the values of the x parameter evolve essen- 
tially identically until ^250 ms and stay below the criti- 
cal threshold of 3 for convection for most of the time. 
Correspondingly, both the 2D and 3D simulations do 
not develop any significant level of convective activity. 
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Therefore pure SASI activity can be observed, which in 
both cases grows in a very similar way and is character- 
ized by a dipolar shock deformation mode that amplifies 
most rapidly and that dominates the I = 2, 3 spherical 
harmonics components in 3D even more than in 2D. 

In model S25-3, which differs from S25-2 only in the 
size of the initial seed perturbations, the quick appear- 
ance of buoyancy asymmetries leads to a stronger shock 
expansion, which prevents x horn dropping below 3 
in both the 2D and 3D runs. In this case SASI ac- 
tivity cannot be unambiguously diagnosed by persis- 
tent shock sloshing or spiraling. The amplitudes of the 
£ = 1, 2, 3 shock deformation modes in 3D remain fairly 
small (much smaller than the maximum amplitudes of 
S25-1 and S25-2), and a clear hierarchy of these lowest- 
order multipoles over longer periods of time is not estab- 
lished in the 3D case of S25-3 (Figure [9j bottom right 
panel) . 

Models S25-2 and S25-3 in comparison confirm that the 
growth conditions of SASI activity improv e in phases of 
strong shock contraction. As discussed in IS check et al.l 
(2008), the SASI growth rate roughly scales inversely 
with the advection timescale of the accretion flow from 
the shock to the neutron-star surface. In model S25-3 the 
stronger shock expansion, triggered by convection that is 
seeded by the (nonlinear) initial perturbations, leads to 
longer advection times and thus disfavors the develop- 
ment of the SASI. 

Model S25-1 demonstrates another interesting effect. 
The SASI can be stronger in 3D models than in the cor- 
responding 2D cases. In the 3D version of S25-1 the shock 
deformation mode develops a very large t = 1 amplitude, 
for which after 250 ms cf exceeds the amplitudes of the 
higher-order multipoles by more than an order of magni- 
tude (Figure El bottom left panel). While quasiperiodic 
shock oscillations with nice regularity and nearly stable 
axis characterize the evolution until about 350 ms, the di- 
rection of the dipole axis wanders afterward to give way 
to a strengthening spiral mode after *~500 ms, which co- 
incides with the phase of most extreme shock recession. 
In the 3D case the SASI activity appears much more 
strongly than in 2D, where dipolar and quadrupolar de- 
formation modes have similar size and higher multipoles, 
associated with convective mass motions, contribute sig- 
nificantly. The relevance of convection in the 2D model 
is compatible with the x being around 3 in this case, 
whereas x < 3 in the 3D run reduces the influence of 
buoyancy but signals more favorable conditions for the 
SASI because of a shorter advection timescale (similar to 
the situation in S25-2). 

It is important to note that both S25-1 and S25-3 were 
computed with large initial seed perturbations. In the 
case of S25-3, however, the larger boundary and accre- 
tion luminosities lead to stronger neutrino heating, which 
supports buoyancy and allows the shock radius to stay 
large even in the 3D case. On the contrary, in S25-1 the 
weaker neutrino heating does not prevent shock retrac- 
tion. 

Convection and SASI activity in 3D therefore depend 
sensitively on the behavior of the accretion shock, which 
in turn reacts to the neutron star contraction and the 
power of neutrino-energy deposition. Our set of mod- 
els shows that for a favorable combination of conditions 
the SASI in 3D can be much stronger (S25-1) and purer 



(S25-2) than in the corresponding 2D cases. In gen- 
eral, however, the development of convection and SASI 
in the postshock accretion layer seems to be fully com- 
patible with the general theoretical understanding ob- 
tained in connectio n of 2D hydrodynamic simulations in 
iScheck et all (pOOl . 

4. DISCUSSION AND CONCLUSIONS 

We have simulated the post-boun c e evo lution of the 
27 M Q progenitor of iWooslev et al.l (|2002| ) in 2D and 
3D, using the Prometheus- Vertex code with detailed 
multi-group neutrino transport including the full, sophis- 
ticated set of neutrino reactions applied previously in 
ID and 2D supernova modeling by the Garching group. 
Moreover, we performed a set of 2D and 3 D post-bounce 
simula tions of the 25 M Q progenitor of IWooslev et al.1 
(2002) with the Prometheus hydrodynamics scheme 
and a computationally efficient, gray neutrino trans- 
port approximation, employing a parametric approach in 
which an inner grid boundary replaced the excised high- 
density core of the proto-neutron star and allowed us to 
more systematically explore the influence of different core 
neutrino luminosities and of faster or slower contraction 
of the forming remnant. We emphasize that this ap- 
proach is different from that for our 27M Q runs, where 
the whole neutron star to the center was included in the 
computational domain. In addition, we also tested the 
effects of varied amplitudes of random perturbations that 
had to be imposed for seeding the growth of nonradial 
hydrodynamic instabilities in the neutrino-heated accre- 
tion layer behind the stalled supernova shock. While the 
3D run for the 27 M star was conducted with a polar 
coordinate grid, an axis-free Yin- Yang grid was used in 
the studies of the 25 M Q progenitor. 

Our simulations for the first time provide unambigu- 
ous evidence of the occurrence of large-amplitude SASI 
shock sloshing and spiral motions and of their interplay 
with neutrino-driven convection in 3D supernova core 
environments modeled with a "realistic" treatment of 
neutrino transport and the corresponding heating and 
cooling effects. Previous hydrodynamic studies had been 
able to identify SASI activity, and in particular SASI spi- 
ral modes, onl y in 3D setups with adiabatic postshock 
accre tion flows (jBlondin fc Mezzacap pa 2007; Fernandez! 
2010) and in some simulations with a simple neutrino 
light-bulb treatment (jlwakami et al.l I2008L 12001 ) . The 
shallow water analogue of the SASI was also observed 
experimentally, ho wever again without the presence of 
buoyancy motions (|Foglizzo et al.H2012D . On the basis of 
more recent 3D simu lations with a simple neutrino light- 
bulb approximation dBurrows et al"i r2012; Mu rphy et al.l 
l2012tlDolence et alJl2013[) and 3D GR m odels with a neu- 
trino leakage scheme (|Ott et akl 120121 ) it was even con- 
cluded that the SASI is at most of minor relevance for 
the dynamics of the postshoc k layer in collapsing stellar 
cores (see also lBurrows][2013l) . 

Our findings indicate that the SASI is potentially much 
more important in 3D than suggested by these previous 
investigations. Besides the core-density profile of the pro- 
genitor star, which determines the temporal evolution of 
the mass infall rate, the more elaborate neutrino trans- 
port treatment (in particular in our 27 Mq simulation) 
may partly be responsible for the different accretion-flow 
dynamics seen in our models. The details of the neutrino 
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description are relevant, because they can cause differ- 
ences in the contraction behavior of the proto-neutron 
star and in the neutrino-heating in the gain layer. Both 
affect the evolution of the stagnation radius of the accre- 
tion shock and thus have a strong impact on the growth 
conditions for con yection and th e SAS I. In the artifi- 
cial se t up used bylBurrows et al.l (|2012t ); iMurphv et al.l 
( 2012); iDolence et al.l (|2013[ ). the neutron star is not al- 
lowed to deleptonize and radiate energy, therefore its 
radius stays unrealistically large (50-60 km) during the 
whole simulated post bounce evo l ution (up to ~1 s). Also 
the leakage scheme of lOtt et al.l <|2012f ) . which can track 
the contraction of the proto-neutron star and the time 
evolution of the neutrino emission to a certain extent, de- 
viates significantly from full transpo rt models. For exam- 
ple, in the work of lOtt et all (|2012f ) the mean energy of 
the radiated electron neutrinos around 100 ms after core 
bounce is more than 60% hig her than that found with 
the Vertex transport code bv lMiiller et al.l (|2012al) . al- 
though both groups investigated the same 27 Mq progen- 
itor with relativistic methods. It is cl ear that the neu- 
trino heating found bv lOtt et al.l (|2012f l was stronger and 
thus more favorable for larger shock radii, which is turn 
was supportive of neutrino-driven convection instead of 
the SASI. 

In contrast t o the fin dings by iBurrows et al.l (|2012t ); 
Murphy et all (120121) ; iDolence et al.l (|2013l ) and by 
Ott et al.1 (|2012h . we observe strong SASI activity in our 
3D simulations both for the 27 Mq progenitor and the 
25 Mq star. The SASI can become the clearly strongest 
nonradial instability during at least some parts of the 
evolution and can be clearly identified by its oscillatory 
growth properties and even in the nonlinear regime by 
its characteristic low-order spherical harmonics modes of 
coherent shock motion and shock deformation. In detail, 
we can draw the following conclusions from our models: 

• SASI activity can develop in 3D despite and in 
addition to the earlier presence of neutrino-driven 
buoyancy. Different from the higher-order multi- 
pole pattern that is typical of convective plumes 
and downdrafts, the SASI reveals itself by coherent, 
large-amplitude shock sloshing and spiral motions. 

• Because of the absence of a flow-constraining sym- 
metry axis, which also directs the structure of 
neutrino-driven buoyancy in 2D models, SASI 
shock motions and deformation modes can be rec- 
ognized more easily and more clearly in the 3D 
case. Interestingly, both the 27 Mq and 25 Mq 
models exhibit evolution phases in which the SASI 
in 3D can become stronger than in the correspond- 
ing 2D runs. The dominance of the SASI and 
greater strength in 3D can be concluded not only 
from the large dipole and quadrupole amplitudes 
of the shock deformation, but also from the higher 
nonradial kinetic energy of the postshock flow (Fig- 
ure [2]) and from a prominent peak at low-order 
multipoles in the power spectrum of the nonra- 
dial kinetic en ergy (Figure [3 . While some authors 
hypothesized (|Iwakami et all 120081) that the SASI 
amplitudes in 3D remain smaller than those in 2D 
because the kinetic energy of the nonradial flow is 
shared with an additional degree of freedom, our 



results suggest that there is no such limitation of 
the kinetic energy that can be stored in lateral and 
azimuthal mass motions. The fraction of the huge 
reservoir of accretion energy that is converted to 
nonradial flows in the postshock layer can be larger 
in 3D than in 2D. 

• Depending on the conditions in the postshock flow, 
which in our 25 Mq runs could be controlled by the 
choices of the contraction behavior and neutrino 
luminosities assumed at the inner grid boundary, 
phases of essentially pure SASI activity could be 
obtained, with neutrino-driven convection only de- 
veloping as the secondary instability (cf. Figure [T0l 
middle panels) . This result of the 2D and 3D mod- 
els of the present work is very similar to the flow 
behavior seen in the 2D simulatio n of the 27 Mq 
progenitor bv lMiiller etaLl (|2012al) . 

• The growth of the SASI is favored by fast advec- 
tion flows, because its growth rate in an advective- 
acoustic cycle scales roughly inversely with the ad- 
vection timescale of the accretion fl ow from the 
shock to the neutron star surface fe.g. JScheck et ahl 
2008). This is opposite to neutrino-driven buoy- 
ancy, whose growth in the accretion flow requires 
that the critical ratio x °f advection timescale to 
buoyancy timescale exceeds a value of about 3 
(Fogli zzo et al.l 120061 ) . Fully consistent with this 
dimension-independent theoretical understanding, 
which was developed by linear analysis and hydro- 
dynamical modeling in 2D, we find SASI growth 
also in the 3D case to be strongest in phases of 
small shock stagnation radii, which are connected 
to rapidly shrinking and more compact proto- 
neutron stars as well as relatively weak neutrino 
heating. 

• Also the amplitude of the initial seed perturbations 
(of the density or velocity field) is found to poten- 
tially have an influence. A greater amplitude can 
trigger a faster development of neutrino-driven con- 
vection (for perturbations of nonlinear size this can 
happen even whenx < 3, see lScheck et~aT1l2008l ). If 
the buoyancy instability is supported by sufficiently 
strong neutrino heating to instigate an expansion 
of the shock, SASI growth can become disfavored 
and nonradial flows in the postshock layer are dom- 
inated by buoyant plumes and downflows. 

• Preferentially in phases of strongest shock recession 
we observe bipolar (£ = 1, m = 0) SASI sloshing 
motions to change over to a spiral (£ = 1, m = 1) 
mode in the 3D simulations of both the 27 Mq and 
25 Mq models (Figures [6] and [7]) ■ The transition 
to the time-dependent rotating shock-deformation 
pattern therefore seems to be favored by particu- 
larly small shock radii. In general, the character 
of the 3D accretion flow during shock oscillation 
and spiraling phases with wandering directions re- 
veals close similarity to the behavior of the hy- 
draulic jump observ ed in the SWASI experiment 
(jFoglizzo et al.H2012f ). 

Our comparison of the SASI activity in the 2D and 
3D simulations of the 27 Mq model, which was evolved 
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fully self-consistently (i.e., without the use of free pa- 
rameters and without an inner grid boundary) is par- 
ticularly interesting. It revealed that at an early stage, 
when the mass accretion rate is still high and the shock 
correspondingly retreats in response to the proto-neutron 
star contraction, the SASI can grow despite some prior 
neutrino-driven convective activity. Before the accretion 
of the Si/SiO interface, the SASI can reach even higher 
amplitudes in 3D than in 2D. The large-amplitude shock- 
sloshing mode eventually makes the transition to a spiral 
mode in 3D, which, however, is quenched around 260 ms 
after bounce once the accretion rate has dropped signif- 
icantly and the shock has expanded to radii of nearly 
200 km on average. Despite the faster shock expansion 
and higher kinetic energy of the postshock flow, the 3D 
model nevertheless falls short of an explosive runaway 
unlike the 2D simulation. 

A discussion of this important difference in the explo- 
sion behavior of 2D and 3D model s, which is fu l ly con - 
sistent with the results obtained bv lHanke et al.l (2012), 
is beyond the scope of this paper and will be addressed 
in future work. We do not want to speculate here about 
the relevance of the SASI in competition with neutrino- 
driven convection for getting explosions in 3D. It is pos- 
sible that the SASI could provide crucial support for the 
onset of the explosion, especially since a very strong spi- 
ral mode returns towards the end of our 27M Q simula- 
tion. Such a SASI supported explosion was obtained in 
our parametric 25 M© simulation of model S25-2 (Fig- 
ure 121 middle column) . But it is also possible that its 
presence is only temporary and that neutrino-driven con- 
vection dominates when the explosion sets in and the pre- 
ceding shock expansion leads to disfavorable conditions 
for the advective-acoustic cycle feeding the SASI. Even in 
the latter case, the SASI could still play an essential role 
in the explosion mechanism by pushing out the shock 
far enough for convection to take over in cases where 
convection cannot develop on its own in an initially sta- 
bilized post-shock flow. The conditions for the growth 
and long-time persistence of the SASI, in particular of 
the spiral mode, are still poorly understood and have 
to be explored in more detail by future 3D simulations. 
Our results suggest that the exact role of the SASI will 
depend strongly and in a subtle way not only on a vari- 
ety of progenitor dependent conditions like the steepness 
of the density profile, the location and sharpness of the 
composition shell interfaces, the existence of initial non- 
radial asymmetries in the collapsing stellar core, angular 
momentum of the progenitor core, which could trigger a 
faster growth of the spiral mode e ven for small rotation 
rates (Yamasaki & Foglizzo 2008), and magnetic fields, 
if these are amplified to dynamically relevant strength. 
The SASI growth conditions will also depend on the 
contraction behavior of the nascent neutron star, which 
is influenced by its high-density equation of state and 
by general relativistic gravity, and on the strength of 
the neutrino heating, whose determination requires a so- 
phisticated treatment of the energy-dependent neutrino 
transport. 

The initial perturbations constitute a potential cause 
for differences that is very difficult to address. Cer- 
tainly neither the artificially imposed random perturba- 
tions i n our models, nor th e small intrinsic numerical 
seeds of lMiiller et a l. (2012a), nor the larger intrinsic per- 



turbations in the simulations of lOtt et al.l (|2012f) reflect 
the magnitude and spectral shape of the physical per- 
turbations present in realistic progenitors. In the end, 
a well-founded understanding of the interplay between 
convection and the SASI will probably require multi- 
dimensional progenitor models. 

If we suppose that the growth conditions for the SASI 
remain favorable for realistic seed perturbations, our self- 
consistent 3D model of the 27 Mq progenitor still throws 
up a number of further questions: Has the SASI already 
reached its saturation level in our model, or could it 
grow further if massive accretion were to last longer? 
Could it play a decisive dynamical role for reaching an 
explosive runaway in the latter case, even if it is only 
by pushing the shock out far enough for convection to 
take over later? Are the saturation properties of the 
nonlinear spiral mode different from the nonlinear slosh- 
ing mode because of the different geometry of the down- 
flows and the presence of considerable angular momen- 
tum in the gain region? If the strong and clear spiral 
mode survives until the onset of an explosion, what would 
be its final impact on the proto-neutron star kick and 
spin? What are the neutrino-emission variations con- 
nected to large-amplitude shock sloshing and spiral mo- 
tions in the 3D case? Are these signal variations observ- 
able wfilijajjnm^ 2D case 
Csee lMarek fc Jankal[2009t ILund et alJlMlCh ? With so- 
phisticated 3D supernova modeling only now beginning 
to become feasible, our understanding of 3D effects in su- 
pernova cores is still in its infancy, and further research 
into the properties of the SASI is urgently needed. 
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